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Abstract 

We investigate nonfactorizable contributions to charm 
meson decays in D — > Kir/ Kp /K*tc / Ka\ / K* p modes. 
Obtaining the contributions from spectator-quark diagrams 
for N c = 3, we determine nonfactorizable isospin 1/2 and 
3/2 amplitudes required to explain the data for these modes. 
We observe that ratio of these amplitudes seem to follow 
a universal value. 
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1 Introduction 



As precise data [1] on some of the weak hadronic and semileptonic 
decays of charm mesons is available, it has now become possible 
to test the validity of the factorization model. Employing isospin 
formalism, in a strong interaction -phase independent manner, re- 
cently, Kamal and Pham [2] has, shown that the naive factorization 
model fails to account for isospin amplitudes ( A 1 / 2 and A 3 / 2 ) for 
D — > Kit, D — > Kp, and D — > K*ir modes. It was observed that 
for D — > PP (P = pseudoscalar meson), whereas the factorization 
assumption accounts for branching ratio of D + — > K n + , it overes- 
timates A1/2 for D° -decays. For D — > PV ( V = vector mesons) 
modes, it overestimates A 3 / 2 for D — > Kp, decays and underesti- 
mates A3/2 for D — > K*ix decays. One of the ways to remove the 
discrepancy could be the inelastic final state interaction [2,3] which 
can feed K*ir channel at the expense of K* p channel. An alter- 
native may be the nonfactorizable contributions [4] arising through 
the soft gluon exchange, which are generally ignored in the factor- 
ization model. Recently, there has been a growing interest [5,6] in 
exploring such contributions in the hadronic decays of charmed and 
bottom mesons. Major cause of the concern has been that N c — > 00 
limit, which is considered to be supported by D-meson phenomenol- 
ogy, fails when carried over to B-meson decays [7]. Therefore, a 
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reinvestigation of the charm decays is called for. 

We, in this paper, study the nonfactorizable contributions to 
D — > Kir, D — > Kp, D — > 7?*tt, D -> X Ql and £> -> i?*p de- 
cays. Employing the isospin formalism in the phase-independent 
manner, we determine these contributions in respective 1/2 and 3/2 
-isospin amplitudes from experiment and search for a systematics in 
the amplitudes. 

In section 2, weak Hamiltonian is discussed. In next sections, we 
study these decay modes one by one. Summary and discussion are 
given in the last section. 

2 Weak Hamiltonian 

The effective weak Hamiltonian for Cabibbo-angle-enhanced decays 
of the charm hadrons is given by 

H w = G F [ci(ud)(sc) + c 2 (sd)(uc)], (1) 

where G F = ^V ud V* s and qiq 2 represents color singlet V - A current 

q iq2 = gi7 M (l - 75)^2, 

and the QCD coefficients at the charm mass scale are [8] 

ci = 1.26 ±0.04, c 2 = -0.51 ±0.05. (2) 
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Due to the Fierz transformation of the product of two Dirac currents 
in (1) in 7V c -color space, the Hamiltonian takes the following form 

H C J = G F [ ai (ud)(sc) + c 2 H 8 w }, 

H° s = G F [a 2 (sd)(uc) + Cl H*}, (3) 

for color favored (CF) and color suppressed (CS) decay amplitudes, 
and 

ai,2 = c 1>2 + — , 
Hi = i^( M A^)( S -A a c), 

Z a=l 

Hi = ^(sA^AV), (4) 

a=l 

where color-octet currents 

q 1 \ a q 2 = qil^l - 75)A a g 2 

involve the Gell-Mann matrices, A a for color. Matrix element of the 
first terms of (3) can be calculated using the factorization scheme 
[9]. In this scheme, nonfactorizable contributions from the second 
terms in (3) are ignored, so long one restricts to color singlet in- 
termediate states. Due to the neglect of these terms, one usually 
treats ai and a 2 as input parameters instead of using iV c = 3 in 
reality. Empirically D — > Kit data seem to favor 7V C — > oo limit 
[9] which is justified with the hope that the nonperturbative effects 
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arising due to the soft gluon exchange between the colored-octet 
current in (4), relative to the factorizable amplitudes arising from 
the color singlet current in (3), are of the order in the large N c - 
limit. In fact, in some of the QCD calculations [10], it has been 
claimed that for B — > Dn, and D — > Kix as well, nonfactorizable 
terms tend to cancel the contributions from jj- terms in the first 
terms of (3). However, B-mesons don't favor this result empirically. 
Further, this does not guarantee that such cancellations would per- 
sist for other decay modes too. In other words, a\ and a 2 are not 
universal parameters, these are decay dependent if one is to stick to 
the factorization model. An alternative is to take N c = 3 and inves- 
tigate nonfactorizable contributions more seriously. This is what we 
do in the following sections employing the flavor-isospin framework. 

3 D — > Kjt decays 

In terms of the isospin amplitudes and final state interactions phases 
[2], 

A(D° - K-n + ) = -L[A 3/2 e^/ 2 + V2A 1/2 e iS ^\, 

A(D° -> K\°) = -L[V2A 3/2 e i5 ^ - A 1/2 e iS ^}, 

v3 

A{D + ^K\ + ) ----- V3A 3/2 e iS ^. (5) 
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These lead to the following phase independent quantities: 

\A(D° - K--k + )\ 2 + \A(D° - K\°)\ 2 = \A 1/2 \ 2 + \A 3/2 \ 2 , 

\A{D+^K\+)\ 2 = 3\A 3/2 \ 2 . (6) 
With Experimental [1] values: 

B(D°^K-7i + ) = (4.01 ±0.14)%, 

B (D° -> XV) = (2.05 ± 0.26)%, 

B(D + ^K\ + ) = (2.74 ±0.29)%, 
and D-meson lifetimes 

r D o = 0.415 ps, 

r D + = 1.057ps, 

decay rate formula 



T(D^PP) = \Gf\ — t-—\A(D — »■ PP 

8nm z 



2 P _| 4/ n . nnM2 



yields 

|^i/2 |«rp = 0.387 ± 0.011 GeV 3 , 

\A 3/2 \ exp 0.097 ± 0.005 GeV 3 (7) 

Since the relations (6) are independent of the final state interac- 
tion phases, one might well evaluate them without the phases and 



determine nonfactorizable contributions using experimental values. 
We separate the factorizable and nonfactorizable parts of the decay 
amplitude as 

A(D -> Kir) = A f (D -> Kir) + A nf (D -> Xtt). (8) 

Using the factorization scheme, factorizable part of the decay am- 
plitudes can be written [9] as 

Af(D°^K-n + ) = aMml-m 2 K )F° K (ml), 

= 0.351 GeV 3 , 
Af(D°^K\°) = -^=a 2 f K (m 2 D -ml)F^(m 2 K ), 
= -0.030 GeV 3 , 
Af{D+^K\+) = aMm 2 D -m 2 K )F DK (ml) 
+a 2 f K (m 2 D -ml)F Q D ^m 2 K ), 

= 0.309 GeV 3 (9) 
Numerical input for these terms is taken as 

ax = 1.09, a 2 = -0.09, 

f w = 0.132 GeV, f K = 0.161 GeV, 

and 

F DX (0) = 0.76, F^(0) = 0.83 (10) 
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from Ref [11]. Using isospin C. G. Coefficients, nonfactorizable part 
of the decay amplitudes can be expressed as scattering amplitudes 
for spurion + D — > K + n process: 

A nf (D° - K-7I+) = l -c 2 {< Kir\\Hl\\D > 3/2 +2 < Kit\\Hl\\D > 1/2 ), 
A"/(£>° - #V) = ^ Cl (< Kix\\Hl\\D > 3/2 - < ^7r||^||L> > 1/2 ), 

A n '(D+ -> ^°7T+) = C 2 < X7r||^||L> > 3/2 +Ci < K>K\\Hl\\D > 3/2 

(11) 

In order to reduce the number of unknown reduced amplitude, we 
assume the following 

< Kir\\H 8 w \\D > 1/2 < K7r\\H 8 J\D > 1/2 , (12) 

< Kn\\H s w \\D > 3/2 = < Kir\\H*\\D > 3/2 , (13) 

as both and behave like |1, 1 > component of an isovector 
spurion. In fact, and transform into each other under the 
interchange of u and s quarks. Thus, in the limit of flavor SU(3) 
symmetry, the constraints given in eqs. (12) and (13) become reli- 
able. From (5), we can write 

A\UD -> Kn) = -^={V2A n f (D° — > K~n + )—A n f (D° — > K n )} 
v3 

A^Ud^Ktt) = -L{A nf (D° ^ K-n + ) + V2A nf (D° ^ R\ )} 
v3 

= .L{A n f(D + -> K\ + )} (14) 
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Relations (12) and (13) then lead to the following prediction: 

Atf 2 (D -> Kn) = cl + 2cj 
A n 3 J 2 (D^Kn) V2(4-cj) 

= -1.123 ±0.112. (15) 

Experimentally, the nonfactorized isospin amplitudes is determined 

as 

A\ s j2 = +0.082 ± 0.032 GeV 3 , 

A n J /2 = -0.081 ± 0.023 GeV 3 , (16) 

from A 1/2 = 0.387 ± 0.011 GeV 3 , and A 3/2 = 0.097 ± 0.005 GeV 3 
given in (7). This yields 

A nf 

= -1.011 ±0.250, (17) 

^3/2 

in good agreement with theoretical expectation (15). Such isospin 
formalism can easily be extended to D — > Kp, D — > K*n, and 
D — > Kdi decays. Since the isospin structure of these decay modes 
is exactly the same as given in (5) and (11), the same value — 1.123± 
0.112 would follow for the respective ratio of the nonfactorizable 
isospin amplitudes A™J 2 and Ay 2 in each of these cases. In the 
following, we determine these amplitudes from experimental values 
of the branching ratios of these decay modes, and compare their 
ratio with theoretically expected one. 
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4 D — > decays 

We begin with the definition of the decay amplitude A(D — > PV) 
through the decay rate formula, 

T(D^PV) = \G F \ 2 ^-\A{D ^ PV)\ 2 . 

A(D — > PV) has the dimension of GeV and is obtained by writing 
down the standard decay amplitude and removing from it a factor 
2m v (e* .P D )Gp, where m v is the vector meson mass and e* is its po- 
larization vector, Pjj is the D-meson four momentum. Factorizable 
parts of the decay amplitudes A(D — > Kp) can be written as 

Af(D°^K-p + ) = ai f p FP K (m 2 p ), 

= 0.203 GeV, 

Af(D°^K°p ) = ^=a 2 f K A^(m 2 K ), 

= -0.0074 GeV, 

A f (D + -> K°p + ) = a,f p F° K {m 2 p ) + a 2 f K A^(m 2 K ), 

= 0.192 GeV. 

Here, we use 

f p = 0.212 GeV, 

F? K (0) = F DK (0) = 0.76, 
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from (10) and 

Aq P (0) = 0.669 

is taken from the BSW model [9]. Numerical values given above are 
calculated by extrapolating F® K (q 2 ) and AQ p (q 2 ) using a monopole 
formulae with pole mass 2.11GeV (D* pole) and 1.865 GeV ( D- 
pole) respectively. Experimental [1] values of branching ratios: 

B(D°^K-p + ) = (10.4 ±1.3)%, 

B(D°^K°p ) = (1.10 ±0.18)%, 
B(D + ^K°p+) = (6.6 ±2.5)%, 
yield the total isospin-amplitudes: 

\A 1/2 \ex P = 0.235 ± 0.014 GeV, 

\A 3/2 \ exp = 0.067 ± 0.013 GeV (19) 
Writing A n f (D — > Kp) analogues to (11) - (13), inturn leads to 

A"f 2 (D -> Kp) = ±0.065 ± 0.015 GeV, 

A^J 2 (D^Kp) = -0.041 ± 0.013 GeV, (20) 



with positive signs chosen for both A±/ 2 and A 3 / 2 in (19). The ratio 
v -1.481 ±0.582, (21) 



A\) 2 {D -» Kp) 
A n 3 j 2 (D - Kp) 



I"/ 

1 3/I 

is consistent with (15) within errors 
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5 D — > X*7r decays 

Repeating the same procedure used for D — > .ft'p decays, here the 
factorizable amplitudes are given by 

A f(D°^K*-n + ) = ai f n A° K \ml), 

= 0.100 GeF, 
Af(D° -> 7TV) = ±=a 2 f K *F?*{m 2 K *), 
= -0.015 GeV, 
A / (jD + _» x*°tt + ) = oi/^'K) + a 2 f K *F^{m 2 K *), 

= 0.080 GeF. (22) 

Here, we use 

F-^O) = F D "(0) = 0.83, f K * = 0.221 Ge\/, 
and Aq K *(0) is determined from the relation 

< x *(0) = A° K *(0) = -^{( mD + mK *)A? K *(0)-(m D -m K *)A° K * 

(23) 

Semileptonic D — > if*ez/ data [11] yields 

Af^(0) = 0.61 ±0.05, A% K '(0) = 0.45 ±0.09. 
Which gives 

A° K *{0) = 0.70 ±0.09. 
12 



Experimental [1] values of branching ratios: 

B(D° -> K*~ir + ) = (4.9 ±0.6)%, 

B(D°^K*\°) = (3.0 ±0.4)%, 
B(D + -> K*°n + ) = (2.2 ± 0.4)%, (24) 

yield the total isospin-amplitude: 

\A 1/2 \ exp = 0.186 ± 0.009 GeV, 

\A 3/2 \ exp = -0.036 ± 0.003 GeV (25) 

Choosing positive and negative sign for the A±/ 2 and A 3 / 2 terms 
respectively, we find 

Ay 2 (D -> K*n) = ±0.096 ± 0.009 GeV, 

A n J /2 (D -> TTtt) = -0.082 ± 0.008 GeV, (26) 

leading to 

11/21 - -1.171 ±0.158, (27) 



A n 3 J 2 (D -> K*tt) 



for monopole like extrapolation of Ff >7T (q 2 ). This ratio is consistent 



with theoretical expectation (15) within errors. Also note that 

\A n J 2 {D - X*tt)| > \A n J 2 {D - Xp)|, (28) 

consistent with theoretical expectations that final state having low 
momentum is likely to be affected more by the soft gluon exchange 
effects. 
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6 D — > Xai decays 

D — > .ft'ai decays can also be treated in a manner similar to that 
used for D — > modes; due to the similarity in their Lorentz 
structure. For D — > Kai decays, factorized amplitudes are ( upto 
the scale factor ^§(e* .p)2m ai ): 

Af(D°^K-ai) = aj^imlj, 

= 0.285 GeV, 
Af{D°^K°al) = -^a 2 f K V D ^(m 2 K ), 
= 0, 

Af{D+ -> ^V) = + a 2 /^ Dai «), 

= 0.285 GeV. (29) 

Here, we use 

/ ai = 0.221 GeV, 

and take 

V DK (0) w 

due to the orthogonality of the D and ai spin-wave functions. The 
experimental values [1,12] for branching ratios: 

B(D° -> X-a+) = (7.9 ± 1.2)%, 

5(D° -> K°a\) = (0.43 ± 0.99)%, 
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B(D + — > K°af) = (8.1 ±1.7)%, (30) 

when used for analogues of relations (6), yield total isospin ampli- 
tudes: 

IA/2U = (0.582^55)^, 

\A 3/2 \ exp = (0.338^;° 6 7 D GeV. (31) 

Here we have neglected the small effects arising due to the width 
of di meson [13]. Choosing positive and negative signs for Ay 2 and 
A 3 / 2 respectively in (31), we find 

A"j 2 (D -> Ken) = +0.349 ± 0.060 GeV, 

A n J /2 (D -> Km) = -0.221 ± 0.023 GeV, (32) 

leading to 



A"L(D -> Xai) 

^_ = -0.910 ±0.165, (33) 

which is consistent with theoretical expectation (15). Also note that 

\Atf 2 {D - PA)\ > \A r ;J 2 (D -> PV)|, (34) 

again in accordance with theoretical expectation, as in the D — > Xai 
decays, final state momentum is smaller than that in D — > PV^ 
mode. 

Thus, we observe that in all the decay modes, considered so far, 
D — > A'tt, P> — > Kp, P> — > X*7r, D — > Kai, the nonfactorizable 
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isospin amplitude Ay 2 not only have the same sign for these decays, 
but also bears the same ratio, with in the experimental errors, with 
A^L amplitude, i.e., 

A\%{D -> Raj) ^ A^D^K**) ^ A\] 2 {D^Kp) 



A»UD - K 



1/ 



7T 



(35) 



Further, we notice that the nonfactorized amplitudes show an in- 
creasing pattern with decreasing momenta available to the final state 
particles, i.e., 

\A nf {D ->#ai)| > \A nf {D ^ K*tt)\ > \A nf {D^Kp)\ 

(36) 

This behaviour is understandable, since low momentum states are 
likely to be affected more through the exchange of soft gluons and 
can acquire larger nonfactorizable contributions. If one takes value 
of the ratio of A™J 2 and Ay 2 in (35) to be -1.123 as obtained in (15), 
and determine A^L from D + decay, one can predict the sum of the 
branching ratios of D° -meson decays in the corresponding mode. 
Following this procedure, we predict 

B(D° -> K-7i + ) + B(D° -> K\°) = 6.30 ± 0.67% 

= (6.06 ±0.30% Expt), 
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B(D° -> K-p + ) + 5(_D° -> i?°p°) = 10-17 ± 3.85% 

= (11.50 ± 1.31% £xpt), 

£?(_D° -> + 5(£>° -> i?*°7T ) = 6.29 ± 1.20% 

= (7.9 ± 2.2%, Expt) 

B(D° -> /Tat) + -B(-D° -> X°a°) = 10.67 ± 2.24% 

= (8.33 ±1.56% Expt). (37) 

All theoretical values match well with experiment. Infact, these 
relations can be expressed in a general form as 

B +B - TD °R h | r, | (V2-a)A^-(l + V2a)At c 

B-+ + B 00 - B 0+ [l + {a-\ j j, 

6t d + A 0+ 

(38) 

with a = A™L/A™L, where subscript — h, 00, 0+ denote the charge 
states of strange and nonstrange mesons emitted in each case. A_+ 
and A^ c denote the factorized amplitudes of D° decays. A 0+ is 
obtained from the D + -decay branching ratio B 0+ , 



5, 



o+ 



10+ 



Td+ x (Kinematic factors) 



7 D — > if*p decays 

In general, -D — > W modes involve Lorentz structure for three 
partial waves: S, P, D waves. Therefore, one may expect nonfac- 
torizable contributions to be present in all of them. Experimentally 
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their partial wave structure has been analysed [1,12]. Data indi- 
cates that S-wave is dominant in the D + — > K*°p + decays. For 
D° — > K*°p° mode, D-wave component seems to exist which in- 
terfere destructively with S-wave. P-wave component of these two 
modes is negligible. However, for D° — > K~p + , data is not clean 
enough to separate these partial waves, though P-wave component 
is found to be small here also. Looking at the experimental status, 
we introduce the nonfactorizable term in S-wave, and relate only the 
S-wave decay branching ratios. Then D — > VV decays also share 
the same isospin structure as given in (5) and (11). The factorizable 
decay amplitudes ( in S-wave ) are given by (upto the scale factor 

A f (D° - K*-p+) = ai (m D + m K ,)m p f p A? K *(m 2 p ), 

= 0.330 GeV 3 , 
A f (D° -> K*°p°) = -^=a 2 (m D + m p )f K *m K *A? p (m 2 K *), 

= -0.030 GeV 3 , 
A f( D + ^ K*°p+) = ai (m D + m K *)f p m K *A° K *(m 2 p ) 
+ a 2 (m D + m p )f K *A® p (m 2 K «), 

= 0.289 GeV. (39) 
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Numerical values given above are determined using the form factors 
Af K (0) already given and A? p (0) = 0.775 as determined in the 
BSW model [9]. 

The decay rate formula [14] reduces to 

T (D - VV) = |G F | 2 ^-(2+( m ^ ~ m ^ - m ^ f)\A{D - VV)\ 2 , 
v 87rm 2 D K K 2m Vl m V2 ' n v n ' 

(40) 

for S-wave. Writing 

A(D^VV) = A f (D ^VV)+A nf (D -+W), 
we obtain 

A n 3 j 2 (D^VV) = +0.11418;°£, (41) 

from experimental value of S-wave branching ratio Bs(D + — > K*°p + ) 
(1.7 ± 1.6)%. Extending the apparent universality (35) of the ratio 
to the D — > VV decay modes, 



A\) 2 {D -> VV) 
jqtiD - VV) 



-1.123 ± 0.112. 



We calculate the sum of S-wave branching ratios of D° — > K*~ p + 
and D° -> ^*°p° decays, 

B S (D° — > K*~ p + ) + B S (D° — > K*°p°) = (14.01?;!)%. (42) 

Subtracting the experimentally known value of 

5 5 (£° -> 7T p ) = (3.0 ± 0.6)%, 
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we predict 

B S (D° -> K*-p°) = (11.0l?;jj)%. 

It is interesting to remark this value satisfy the following relation: 

B total (D°^K*-p+) B total (D°^K*°p ) 
B S (D° -> K*-p+) ' B S (D° -> AT*°p°) ' 

0.54 ±0.25 = 0.53 ±0.17, 

indicating the destructive interference of S- and D-wave partial waves 
for D° — > K*~p + decay also. 

8 Summary and Discussion 

In this work, we have investigated the nonfactorizable contributions 
to various decays of D° and D + mesons Kn /K p/K*7r/ Ka\j K* p 
states involving isospin 1/2 and 3/2 final states. In our analysis, 
we take the real value of N c = 3. Since the nonfactorizable con- 
tributions, being nonperturbative, are difficult to be calculated, we 
determine their amplitudes in these isospin states required by the 
experiment. We have ignored the annihilation contributions here as 
these don't contribute to D + decays and for D° -decays these are 
suppressed due to the small value of a2(~ — 0.09±0.05). We observe 
that not only the nonfactorizable isospin amplitudes A™J 2 have the 
sign for the modes considered, but also bears the same ratio with 
^3/2 within experimental errors. The ratio can be understood on 
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the basis of V-spin symmetry, under which nonfactorizable parts 
if® and if® of the weak Hamiltonian transform into each other. 
Further the nonfactorizable contributions also show an increasing 
pattern with decreasing momentum available to the final state par- 
ticles emitted in these decays. Extending the apparent universality 
of the ratio A^L/A^L to D — > VV modes, we predict the S-wave 
branching ratio for B S (D° — > K*~p + ) = (11.0l?;g)%, indicating de- 
structive interference between S-wave and D-wave components for 
this decay. 
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